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ABSTRACT 

We consider cosmological applications of galaxy number density correlations to be inferred from 
future deep and wide multi-band optical surveys. We mostly focus on very large scales as a probe 
of possible features in the primordial power spectrum. We find the proposed survey of the Large 
Synoptic Survey Telescope may be competitive with future all-sky CMB experiments over a broad 
range of scales. On very large scales the inferred power spectrum is robust to photometric redshift 
errors, and, given a sufficient number density of galaxies, to angular variations in dust extinction and 
photometric calibration errors. We also consider other applications, such as constraining dark energy 
with the two CMB-calibrated standard rulers in the matter power spectrum, and controlling the effect 
of photometric redshift errors to facilitate the interpretation of cosmic shear data. We find that deep 
photometric surveys over wide area can provide constraints that are competitive with spectroscopic 
surveys in small volumes. 

Subject headings: cosmology: theory — dark matter — galaxies:clusters:general — gravitational lens- 
ing — large-scale structure of universe 



1. INTRODUCTION 

Results from the Wilkinson Microwave Anisotropy 
Probe (WMAP, Bennett et al. 2003) have revealed a 
CMB temperature power spectrum that is remarkably 
well-fit by a simple 5 or 6-parameter model (e.g. Spcrgcl 
et al. 2003). The success of this simple model is some- 
what qualified by irregularities on large scales. There 
is evidence for departures from statistical isotropy, an 
anomalously low quadrupolc amplitude, an anomalous 
absence of correlations on angular scales greater than 
60 degrees, non-Gaussianity, and sharp features in the 
temperature power spectrum (Luminet et al. 2003; Peiris 
et al. 2003; de Oliveira-Costaet al. 2004; Efstathiou 2004; 
Eriksen et al. 2004b; Hansen et al. 2004; Schwarz et al. 
2004; Land & Magueijo 2005; Jaffe et al. 2005). While 
these irregularities may eventually be revealed to be due 
to improper modeling of the instrument or astrophysical 
foregrounds, they certainly have served to call attention 
to the possibility of interesting departures from the stan- 
dard paradigm on the largest observable scales. In this 
paper we explore how well these scales can be probed 
by forthcoming deep and wide photometric redshift sur- 
veys such as the survey of the Large Synoptic Survey 
Telescope 1 (LSST). 

It has been suggested that by combining CMB obser- 
vations and the Sloan Digital Sky Survey (SDSS) one 
can probe very large scales and thus detect possible fea- 
tures in the primordial power spectrum (Wang, Spergel, 
& Strauss 1999). We find that very large- volume pho- 
tometric redshift surveys can probe these large scales 
with smaller statistical errors. For a sufficiently large 
volume the galaxy power spectra not only provide a com- 
plementary look at structure on very large scales, but 
they can also determine the power spectrum on a given 
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scale even more precisely than with CMB temperature 
anisotropy measurements. The reason is that in deter- 
mining a power spectrum one is fundamentally limited 
by the number of modes on a given scale, and there are 
more such modes with a three-dimensional survey than 
with a two-dimensional survey that encloses the three- 
dimensional survey. 

Determining the power spectrum of the fluctuations 
on very large scales is of great interest for several rea- 
sons. First and foremost, the power spectrum on these 
scales can be cleanly used to infer the power spectrum of 
primordial fluctuations, and the primordial fluctuation 
power spectrum on any scale is of interest. The primor- 
dial power spectrum is one of our only handles on the 
mechanism that led to the fluctuations that are respon- 
sible for the diversity of structures we see in the Universe 
today, including ourselves. 

In the context of inflation, probing larger scales means 
probing inflation at an earlier epoch. Larger scales exited 
the fixed-size horizon earlier and thus had more time 
to expand. Reaching back to these earlier epochs may 
provide us with valuable clues about inflation. 

Finally, exploring these large scales is perhaps the best 
way, other than precision measurement of CMB polar- 
ization on large angular scales, to test the claims of sta- 
tistical anisotropy in the WMAP temperature maps. 

Others have suggested ways to follow up on the large- 
scale WMAP irregularities. To shed further light on the 
largest-scale irregularity of all, the low quadrupole am- 
plitude, Dore, Holder, & Loeb (2004) and Skordis & Silk 
(2004) have proposed to use polarization data. Kesden, 
Kamionkowski, & Cooray (2003) have proposed to use 
all-sky cosmic shear measurements. We must also note 
that the low quadrupole on the sky is actually not that 
unlikely in a ACDM model, and could simply be a sta- 
tistical fluctuation (Bridle et al. 2003; Cline, Crotty, & 
Lcsgourgues 2003; Contaldi et al. 2003; Gaztahaga et al. 
2003; O'Dwyer et al. 2004). More conservative fore- 
ground modeling also reduces the level of discrepancy 
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(Erikscn ct al. 2004a; Slosar & Scljak 2004). There are 
also statistical analyses with results that are consistent 
with statistical isotropy (Hajian, Souradeep, & Cornish 
2005). 

When independent measurements of the cosmic shear 
and galaxy power spectrum are combined, one can ob- 
tain robust constraints on dark energy (Hu & Jain 2004) . 
Furthermore, Pen (2004) showed explicitly that with a 
determination of the galaxy power spectra as a function 
of rcdshift, one can then indirectly infer the cosmic shear 
power spectrum with a lower statistical uncertainty than 
one can directly from the shear maps. While the LSST 
survey will sharply address dark energy via cosmic shear 
tomography of billions of galaxies, here we emphasize 
using galaxy photometry (not shear) in that survey to 
determine the matter power spectrum. 

Deep and wide surveys are not only critical for mea- 
suring the power spectrum on very large scales, they 
also lead to reduced sample variance errors for power 
spectrum measurements on smaller scales. We therefore 
consider applications of intermediate-scale constraints on 
the power spectrum as well. Namely, the broadband 
shape and baryon acoustic oscillations (Peebles & Yu 
1970; Bond & Efstathiou 1984; Holtzman 1989; Hu & 
Sugiyama 1996) in the matter power spectrum can serve 
as standard rulers to determine the angular diameter 
distance r(z) and to constrain cosmological parameters 
(Eisenstein, Hu, & Tegmark 1998; Cooray et al. 2001; 
Hu & Haiman 2003; Linder 2003; Seo & Eisenstein 2003; 
Matsubara 2004). We find the LSST survey capable of 
~ 1% distance measurements to redshifts between 0.2 
and 3. 

The paper is organized as follows. Section 2 describes 
our fiducial survey: the proposed LSST survey. The er- 
rors on the matter power spectrum are forecast in Section 
3, where effects of photometric redshift errors, redshift 
distortion, and spherical survey geometry are illustrated. 
Wc demonstrate in Section 4 that galaxy counts from the 
fiducial survey can be used to determine the differential 
extinction and photometry errors so that these errors do 
not significantly contaminate the power spectrum. Sec- 
tion 5 shows that comparable precision on the matter 
power spectrum can be achieved by the LSST and CMB, 
and that one can measure angular-diameter distances to 
percent level precision from baryon acoustic oscillations 
with the LSST. Wc conclude in Section 6 with remarks on 
several challenges to measuring the matter power spec- 
trum from photometric galaxy rcdshift surveys. A brief 
summary of the spherical harmonic analysis, which is 
convenient for the survey geometry, is given in the Ap- 
pendix. 

2. FIDUCIAL SURVEY 

The LSST will image 23,000 square degrees at high 
galactic latitude deeply in 5-6 wavelength bands from 
0.3-1.1 microns (grizy system). Each band will have sev- 
eral hundred 15-second exposures. Given its wide cov- 
erage and survey depth, the LSST is ideal for measur- 
ing the galaxy power spectrum on very large scales. To 
estimate the galaxy redshift distribution for the 23,000 
square degree survey we degrade the Hubble Deep Field 
North (HDF-N), which has well- measured redshifts, to 
match LSST depth and image quality. 

We approximate LSST's filter system using deep HST 



imaging and ground-based J-band (A cen t ra i ~ 13000 A) 
imaging from the 4-m KPNO telescope. First, we con- 
volve the HDF-N UBVI space images (Williams et al. 
1996) with a 0.8" FWHM Gaussian to simulate the worst 
case seeing conditions. Then, we re-pixelize, add noise, 
and catalog the images to match the expected data qual- 
ity for the final full-depth stack of 500 x 15s coadded 
exposures. For example, the final stack will go to 26.7 
mag and 25.4 mag (10 and 30ct) in the i band. The J- 
band image is left unchanged, because it has the same 
resolution and worse seeing than expected for LSST. 

We then compute photometric redshifts for all de- 
tected objects using a technique based on SED fitting, 
and on a magnitude prior (Margoniner 2005, in prepa- 
ration). The density of objects is ~ 250,000 per square 
degree. Because we need galaxies with well-determined 
redshifts, we choose to keep only the ones with a sharp 
and well-defined peak in its redshift probability distribu- 
tion; i.e., we require that the width of the peak be less 
than 0.04(1 + z). The resulting rms photometric redshift 
error is somewhat larger than this threshold width of the 
rcdshift peak. Our final object density after this cut is 
~ 130, 000 per square degree, with a redshift distribution 
7j g (z) parameterized as: 

n g (z) = 530z 2 e" 2/0 32 arcmin^ 2 . (1) 

This exercise results in an rms of a z = 0.065(1 + z). 
We take this as a conservative upper limit to the rms 
and adopt a z = 0.04(1 + z) as our fiducial value, which 
may be achievable by further color-based cuts or priors, 
although this has yet to be demonstrated. Other surveys 
have achieved similar accuracy. For example, with lumi- 
nous red galaxies from the SDSS, Padmanabhan et al. 
(2005) find that a z ~ 0.03 for z < 0.55 and - 0.06 for 
z < 0.7. We also assume that the photometric rcdshift 
bias is calibratable to 0.01 (1+z) so that it does not sig- 
nificantly degrade the reconstruction of the matter power 
spectrum 2 . This bias error is actually less stringent than 
what is required for precise measurements of dark en- 
ergy equation of state parameters through weak lensing 
or baryon acoustic oscillations. 

3. MEASURING THE MATTER POWER SPECTRUM 

On large scales, the cosmic density field can be approx- 
imated by a Gaussian random field. At redshift z, the 
variance of the Fourier transform of the galaxy number 
density contrast is 

(« 8 (k)5*(k)) = P g (k) = &W(z)P(fc) +n~\ (2) 

where P g (k) is the galaxy power spectrum, P(k) the mat- 
ter power spectrum at z = 0, b(z) the galaxy bias, and 
g(z) the growth factor. We have ignored redshift dis- 
tortion in equation (2) for the moment. If the factors 
multiplying P(fc) in the above equation are known, each 
Fourier mode of the density contrast can be used to make 
a very rough estimate of P(k). 

2 To roughly assess the effect of the redshift bias error, one may 
substitute the bias with a fractional distance error of the same 
magnitude. As shown in Fig. 3, a distance error of 3% (caused 
by an error in the matter fraction O m ) does not degrade the re- 
construction much. Hence, a redshift bias of 1% or less will be 
tolerable. 
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Fig. 1. — Procedures for measuring the large-scale matter power 
spectrum with the LSST survey. 



An unknown scale-dependent and stochastic bias will 
limit our ability to determine the matter power spec- 
trum. However, the galaxy bias is thought to be scale- 
independent and deterministic on large scales. This ex- 
pectation has been confirmed for 0.02 /iMpc" 1 < k < 
0.1/iMpc -1 from the power spectrum analysis of the 
SDSS galaxies (Tegmark et al. 2004), although the bias 
increases with luminosity. Extensive studies with the 
halo model, weak lensing, and simulations (Peacock & 
Smith 2000; Hoekstra et al. 2002; Hu & Jain 2004; Wein- 
berg et al. 2004; Seljak et al. 2005) will help us better 
understand the limits of galaxy bias. Here, we assume 
that the bias is known and scale- independent. 

3.1. Procedures 

Measuring the three-dimensional power spectrum on 
giga-parscc scales brings about two unique issues. 

First, one needs the comoving distance r(z) and co- 
moving angular diameter distance to convert the galaxy 
number density in rcdshift and angular coordinates to 
that in distance coordinates. We assume hereafter that 
the curvature is known so that only one distance is 
needed. For a shallow survey at low rcdshift, Hubblc's 
law suffices to provide an approximation of r(z), but for 
a survey that reaches z = 2.5 one must assume a cos- 
mological model to calculate r(z), or, possibly, obtain it 
from baryon acoustic oscillations. 

Second, because galaxies arc observed on a very long 
light-cone, there is considerable evolution of the bias and 
density fluctuations within the survey. If not accounted 
for, this evolution will contaminate the inferred power 
spectrum on scales (along the light-cone) over which the 
bias or growth factor changes appreciably. 

For a photometric survey, there is yet another issue, 
namely, the photometric redshift error. It suppresses the 
power spectrum in a similar way as the pairwise velocity 
dispersion but with greater magnitude. 

To summarize, we outline in Fig. 1 the procedures of 
measuring the matter power spectrum with the LSST 
survey. The survey provides the galaxy number den- 
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Fig. 2. — Effect of depth (z ma x) on error forecasts for mea- 
surements of the matter power spectrum with the LSST. A 23,000 
square degree photometric galaxy redshift survey is assumed. The 
solid line is the fiducial model power spectrum, while the dotted 
line is the power spectrum generated by the step inflation poten- 
tial (Peiris et al. 2003). The error bars are lcr statistical errors of 
the power spectrum measured in non-overlapping logarithmic bins 
with bin width Afc ~ 0.16fc. The inner error bars are based on the 
simple mode-counting in a cubic volume, e.g. equation (4), while 
the outer ones count spherical harmonic modes using equation (7). 
All the power spectra are scaled to z = 0. 



sity n g (# x ,#y, z), which is converted to n g (x) with the 
knowledge of r(z) from a specific cosmological model de- 
termined by the CMB or from baryon oscillations. The 
raw three-dimensional galaxy power spectrum |-D; n | 2 can 
then be calculated. The subscripts In are the result 
of the spherical harmonic analysis (see the Appendix). 
This raw power spectrum is strongly suppressed in the 
radial direction. With a proper modeling of the effect 
of photometric redshift errors, redshift distortions, and 
galaxy bias, one can estimate cosmological parameters 
(e.g. Cole, Fisher, & Weinberg 1994; Ballinger, Peacock, 
& Heavens 1996; Matsubara & Suto 1996; Matsubara 
& Szalay 2001). After the above correction, one ob- 
tains an estimate of the three-dimensional matter power 
spectrum P/ n . Cosmological parameters can also be esti- 
mated from Pi n through the features imprinted in it and 
be fed back to the beginning of the procedures. Finally, 
the primordial matter power spectrum Pj(fc) can be re- 
constructed from Pi„ with inputs of the matter density 
w m and baryon density toh and compared with that from 
the CMB (e.g. Bridle et al. 2003). 

3.2. Forecasting Errors on P(k) 

For simplicity, and because we are only interested in 
the largest scales, we drop the negligible shot-noise term, 
as justified below. Thus, the standard error in P(k) esti- 
mated from modes within a band of width Afc is simply 



ap(fc) = v /27iVfeP(fc), (3) 

where Nk is the number of independent modes in a shell 
of width Afc centered at fc. For a cubic survey with vol- 
ume V, Nk is given by 

N k = fc 2 AfcV/27r 2 . (4) 
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Variants of these equations can be found in Feldman, 
Kaiser, & Peacock (1994) and Tegmark (1997). We will 
see that this simple description of the errors is a good 
approximation. 

To probe the very largest length scales, we want to go 
as deep as possible. As an example, we show in Fig. 2 the 
power spectrum of a concordance model consistent with 
the WMAP (Spergel et al. 2003) along with the power 
spectrum generated by a step inflation potential (Peiris 
et al. 2003). This latter power spectrum provides a better 
fit to the WMAP data. If we cut the fiducial survey 
at z max = 1-0, then the sample variance in the band 
power is too large to detect or rule out the step potential. 
However, it becomes possible when z max extends to 2.5. 

Our description of the errors on P(k) assumes a cu- 
bic geometry without redshift distortions or photometric 
rcdshift errors. In reality, the survey volumes are wedges 
from a sphere, with rcdshift distortions and errors in the 
radial direction. Spherical harmonic analysis (for details, 
see the Appendix) is well suited for the real situation. A 
simple estimator of the matter power spectrum Pi n can 
be constructed from spherical harmonic modes Di mn of 
the observed three-dimensional galaxy distribution: 



1 = 10 z = 2.5 

max 



Pin = 



1 

Eh 



1 



2Z + 1 ^ 



\ " | n |2 



Ni„ 



(5) 



where I and m enumerate the usual angular modes of 
the spherical harmonics Yj m (f), n is associated with the 
radial modes fc;„ of the spherical Besscl function ji(ki n r), 
and Ni n is the shot noise. The ensemble average of the 
estimator is 

(Pin) = T7~ {(\D lmn \ 2 ) - N ln ) = -Lj2 U L'P(kln'), 



El n 



(6) 

where the ensemble average (\Di mn \ 2 ) is independent of 
m, Uf nn , accounts for photometric rcdshift errors and the 
linear redshift distortion, and Ei n = J2 n ' ^inn' ■ We ne- 
glect the nonlinear redshift distortion because it behaves 
in a similar way as the photometric redshift error but has 
a much smaller effect. 

We calculate the three-dimensional redshift-space 
galaxy power spectrum (|An| 2 ) = (\Di mn \ 2 ) using equa- 
tion (A5) and recover the real-space matter power spec- 
trum (Pin) using the simple estimator equation (6). The 
results are shown in Fig. 3 for Z = 10 and z max = 2.5. 
The redshift-space galaxy power spectrum (dotted lines) 
is strongly suppressed by photometric redshift errors at 
k > H(z)/ca z , 3 while the simple estimator (Pt n ) recov- 
ers the matter power spectrum reasonably well at k < 
0.02 /iMpc -1 assuming an rms photometric redshift error 
of <7 z o = 0.04. The discrepancies between the recovered 
and actual matter power spectra at k > 0.02 h Mpc -1 
are due to the broadening of the window function Uf nn , 
at high wavenumbers and our use of the very simple es- 
timator. More sophisticated estimators can be designed 
to improve the accuracy of the recovery. For <t z q = 0.02, 
there will be no discernible difference between the recov- 
ered and the true matter power spectra within the range 
of Fig. 3. 

3 We drop the subscripts of the wavenumber for convenience. 



10 5 



o 



a. 



10 4 



1 1 1 1 1 

1 

- 


i i i i i i i i 


° ( P ln) - 


- (LDJ 2 ) <7 = 

\ i In i ' z - 


n /P \ wi + Vi 


wrong Q m 




s 


/ 

/ .•■■l-.lj 


\ 


\ 




~J \ 




#\ x 








=\0.02 \* . X " 


0.'Q4 


\* 6 » 


0.07 
i i i ',i i i 1 '• 


\ N 
i- i i i i \ i i 



0.01 
k (hMpc- 1 ) 



0.1 



Fig. 3. — Estimated three-dimensional matter power spectrum 
(Pin) (circles) for / = 10 and z ma x = 2.5. The matter power 
spectrum is estimated for the case cr z o = 0.07, where the rms pho- 
tometric redshift error is parametrized as <r z = ct z o(1 + z). The 
recovery works better with a smaller <t z o. The arrows indicate 
errors of the matter power spectrum reconstructed with a galaxy 
bias that evolves 30% slower than the fiducial bias model. The 
squares show the estimated matter power spectrum under an in- 
correct assumption of the matter fraction Q m = 0.3 instead of the 
fiducial value f2 m = 0.27. The solid line is the real-space matter 
power spectrum that one tries to recover. The dotted lines are 
the redshift-space galaxy power spectrum with the effect of photo- 
metric redshift errors and the linear redshift distortion, while the 
dashed line is the redshift-space galaxy power spectrum with only 
the linear redshift distortion. 



Since the correction is so large, one needs to know 
a z accurately to recover the power spectrum. In the 
(unachievable) case of no rcdshift errors; i.e., a z = 0, 
(\Din\ 2 ) is boosted relative to P(k) by a factor that 
approaches a constant at larger wavenumbers, where 
the plane-parallel approximation is valid and the Kaiser 
(1987) formula applies (Heavens & Taylor 1995). 

Uncertainties of the galaxy bias and cosmological pa- 
rameters can cause systematic errors in the estimated 
matter power spectrum. The errors arise in two ways: 
1) mismatched weight function w(r,k) [for measuring 
the observables, see equation (Al)] and the factor Ei„ 
[for recovering the matter power spectrum from the ob- 
servables, see the estimator equation (6)], and 2) wrong 
length scales in all the measurements and calculations if 
an incorrect cosmological model is assumed. 

Current measurements of the galaxy bias at low red- 
shift have errors of roughly 10% (e.g. Hockstra et al. 
2002; Verde et al. 2002; Seljak et al. 2005). If the bias is 
overestimated or underestimated by 10% at all rcdshift, 
then the normalization of the recovered matter power 
spectrum will be changed by roughly 20% in the oppo- 
site direction without affecting the shape. However, the 
shape can be altered if an error is introduced to the evo- 
lution of the bias. The arrows in Fig. 3 illustrate the 
distortion of the recovered matter power spectrum due 
to a 30% underestimation of the bias growth rate with re- 
spect to the fiducial bias model. This case clearly shows 
that understanding the galaxy bias is crucial to a success- 
ful reconstruction of the matter power spectrum from a 
galaxy survey. 
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Fig. 4. — Similar to Fig. 3, but we show the three-dimensional 
redshift-space galaxy power spectrum (|D; n | 2 ) and recovered mat- 
ter power spectrum (Pi n ) for different multipoles. From left to 
right, the redshift-space galaxy power spectrum (dotted lines) and 
shot noise Ni n (dashed lines) are calculated with I = 1,6,20,40, 
and 70. The matter power spectrum (Pi n ) is estimated with 1 = 1 
(squares) and 70 (circles). The redshift-space galaxy power spec- 
trum increases with the multipole number on scales dominated 
by photometric redshift errors, because at the same wavenumber 
higher multipoles contain a smaller fraction of radial modes, which 
are affected by photometric redshift errors, and a larger fraction of 
angular modes, which are not. 



Hu & Jain (2004) demonstrate that with a 4,000 square 
degree next generation ground based weak lensing survey 
one can achieve percent level constraints on the galaxy 
bias (even as a function of scale) by a joint analysis of 
the shear-shear, galaxy-shear, and galaxy-galaxy corre- 
lations. Hence, we are optimistic that, with the larger 
LSST survey, the error in the galaxy bias will not be the 
dominant source of error in reconstructing the matter 
power spectrum. 

If one assumes a matter fraction Cl m = 0.3 instead of 
the fiducial value fi m = 0.27, then the comoving dis- 
tance to z = 2.5 will be reduced by 3%. This leads to 
roughly a 9% reduction in the raw galaxy power spec- 
trum and an increase of the wavenumber by 3%. We 
recover the fiducial-model matter power spectrum by as- 
suming f2 m = 0.3. The results are shown as squares in 
Fig. 3. The scaling of the matter power spectrum due 
to the change in distances is visible but somewhat com- 
pensated by errors in the factor Ei n and weight function 
w(r,k). The error in f2 m does not produce features in 
the recovered matter power spectrum, so that one can 
still detect primordial features on very large scales. 

Examples of different multipoles are given in Fig. 4. 
One sees that the redshift-space galaxy power spectrum 
(dotted lines) increases with the multipole number on 
scales dominated by photometric redshift errors. The 
reason is that at the same wavenumber higher multipoles 
contain a smaller fraction of radial modes, which are af- 
fected by photometric redshift errors, and a larger frac- 
tion of angular modes, which are not. From the difference 
between multipoles of the redshift-space galaxy power 
spectrum one may actually quantify the effect of photo- 
metric redshift errors (and redshift distortions) without 
knowing <j z . This is analogous to using the quadrupole- 
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Fig. 5. — The window functions Uf i that relate the three- 
dimensional real-space matter power spectrum to the redshift-space 
galaxy power spectrum [see equation (A5)]. They are calculated 
here with £ = 10. Symbols also mark the wavenumbers of several 
discrete modes fc; n ; . They become broader at higher wavenumbers 
because of photometric redshift errors and redshift distortion. The 
window functions of other multipoles are similar to those of I = 10 
but shifted toward higher (lower) wavenumbers for higher (lower) 
multipoles. 



to-monopole ratio to determine the redshift distortion 
parameter and pairwise velocity dispersion (e.g. Cole 
et al. 1994; Peacock et al. 2001). The estimated mat- 
ter power spectrum (Pi n ) with I = 1 fails to recover the 
true power spectrum at k > 0.02 h Mpc -1 because the 
window function U\ nn , is too broad there to guarantee 
the accuracy of the simple estimator (note that Uj Qnn , 
remains narrow). This can be improved with a better 
estimator. 

Fig. 4 also demonstrates that the shot noise Ni n is in- 
deed negligible on scales of interest for the fiducial survey 
with z max = 2.5. Since n g (z) increases rapidly with de- 
creasing redshift, the shot noise for a lower z max is also 
negligible. The rise of the shot noise at small wavenum- 
bers is due to the radial variation of the selection func- 
tion as well as the wavenumber-dependence of the weight 
function [see equation (A6)]. Otherwise, Ni n would be 
wavenumber- independent . 

Given that the shot noise is very low, we can divide 
the galaxies into sub-samples of different luminosity and 
measure the power spectrum of each sub-sample sepa- 
rately. Although sub-sampling cannot reduce the cosmic 
variance, one docs gain by reducing the inhomogencity in 
each sub-sample that is used as a (biased) tracer of the 
cosmic density field. By comparing the power spectra of 
luminosity classes, we can infer the relative bias between 
them (see e.g. Norberg et al. 2001; Zehavi et al. 2002) 
and test the assumption of a scale-independent bias on 
large scales. 

The estimator Pi„ can be further binned to reduce 
the error variance. We calculate the window functions 
Uf nn , that relate the three-dimensional real-space mat- 
ter power spectrum to the redshift-space galaxy power 
spectrum using equation (A6). The I = 10 results are 
shown in Fig. 5. Since the window functions are fairly 
narrow in k space, each P\ n is nearly independent of each 
other. Thus, equation (3) is still valid for forecasting the 
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errors, but the mode counting becomes 



N, 



/sky X)(2Z + l)N l: 



(7) 



i>0 



where / s k y is the fraction of the sky covered by the survey, 
and Ni is the number of radial modes of multipolc I that 
fall in the band. We have neglected the shot noise and 
approximated the reduction in the number of modes due 
to partial sky coverage by an overall reduction of modes 
by a factor of /sky For a real survey, one can incorporate 
the sky cut as well as an angular selection function (or 
mask) numerically. We exclude the monopolc modes, be- 
cause they can be confused with the radial selection func- 
tion. Further modifications to equation (7) are needed 
in the regime where the window function Uf , starts to 
broaden. 

After counting spherical harmonic modes, we find that 
the spherical geometry only mildly increases the sample 
variance of P(k) on the scales of interest, which can be 
identified with the outer error bars in Fig. 2. Hence, the 
simple mode counting with a cubic geometry can still 
serve as a reasonable approximation. 

4. EXTINCTION AND PHOTOMETRY ERRORS 

On very large scales, the variance of density fluctua- 
tions in logarithmic k bins is very small, e.g A 2 (fc) = 
k 3 P(k)/2Tr 2 ~ 1CT 3 at k = 0.01 /iMpc -1 . An unknown 
varying extinction over the wide survey area can cause 
fluctuations in galaxy counts that may swamp the sig- 
nal. If the logarithmic slope of galaxy counts n g (< m) 
as a function of magnitude is s = dlogn g /dm, then the 
fractional error in galaxy counts is 



5 n„ 



In 10 s 5m 



2.5 SJ , 



(8) 



where df // is the fractional error in flux caused by, e.g., 
extinction correction residuals. Observationally s varies 
from 0.6 at blue wavelengths to 0.3 in the red (e.g. Tyson 
1988; Pozzetti et al. 1998; Yasuda et al. 2001), and tends 
to be flatter for fainter galaxies (Metcalfe et al. 2001; 
Liske et al. 2003). To keep this systematic angular fluc- 
tuation well below A(fc), one has to reduce the flux error 
to 1% or better over the whole survey area. This is a very 
conservative estimate, because the power spectrum re- 
ceives contributions from not only angular clustering but 
also radial clustering of galaxies on large scales, which is 
much less affected by the extinction or photometry er- 
rors. 

Extinction A\ (in mag) is related to reddening via 
color excess E(B — V). Aside from a zero-point differ- 
ence of 0.02 mag, the reddening maps made with differ- 
ent methods by Burstein & Heilcs (1982) and Schlcgcl, 
Finkbcincr, & Davis (1998) agree with each other up 
to an rms error of 0.007 mag (Burstein 2003). With 
Ab = 4.3.B(5 — V), the relative error in reddening trans- 
lates to an error of 0.03 mag in B magnitude or a flux 
error of 3%. H I data and galaxy counts were combined to 
produce the reddening map in Burstein & Heiles (1982). 
At that time galaxy counts were typically less than 100 
per square degree and thus prone to statistical error. 
With 250,000 galaxies per square degree 4 , the LSST will 

4 Since the two-dimensional projection of galaxy counts does 
not require redshift information, one can use all galaxies that have 
good photometry. 



achieve the same accuracy in differential extinction with 
two-dimensional galaxy counts alone, as long as the dust, 
or anything that alters the flux, is largely confined in our 
galaxy, or, at least, at very low redshift. 

The number of galaxies N g within an angular window 
9(f) is given by 



(9) 



N g = / n g (r)6(f)d J r 



= J n g (r)e(f)[<5 g (r) + l]d 3 r. 



The mean number counts is 

N g = J n g (r)9(r)d 3 r, 
and the variance is 



a k=Yl J P(k)\ni(k)ei m \ 2 k 2 dk, 



(10) 

(11) 



(12) 
(13) 



where 

e lm = J e(f)y^(f)df. 

We have used the identities 

h g (r)r 2 dr = n g (z)dz 
e ikr = 47r5^i'i,(Ar)y Im (r)y,_ m (k) 

lm 

to derive equation (11). The variance of N g is analogous 
to the rms density fluctuation within a volume (e.g. Zhan 
& Fang 2003). 

For demonstration purpose, we calculate the variance 
with an angular window function 

@{e,4>) = e{6) = e - e ' 2 / 26 ™, 

where 9 and <f> are, respectively, the polar and azimuthal 
angles. The multipolc number I can be approximately 
related to the characteristic scale of the window function 
M by I ~ 360 o /6»m. 

For the effect of the galactic extinction, we use the 
same axial symmetric window function, and place it 
randomly on the reddening map (Schlegel et al. 1998) 
with a restriction that its central galactic latitude |6 C | > 
20° + 1.5f?M- The rms fluctuation of B-band galaxy 
counts within the window function is calculated with the 
conversion Sn g /n g ~ 8Ab = 4:.3SE(B — V). 

The results are shown in Fig. 6. The galactic extinction 
apparently has a dominant effect over galaxy clustering. 
However, the rms fluctuation due to galaxy clustering 
does not include radial clustering because of the projec- 
tion in redshift direction. When the radial information 
is restored, the clustering variance will be much higher. 
If one restricts the window function to higher galactic 
latitude, the rms fluctuation due to extinction will be 
smaller. In addition, one can also correct for the fluctua- 
tions in the extinction; what will hinder the measurement 
of the galaxy power spectrum is rather the uncertainties 
in the extinction. 

If one attributes all the angular variations in the pro- 
jected galaxy surface number density to extinction (or 
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Fig. 6. — Rms fluctuations of galaxy counts due to the extinc- 
tion (dotted line), galaxy clustering (solid line), and Poisson noise 
1/ \J N s (dashed line) within an angular window of size 9y[ . The 
histogram shows the contribution to the number of independent 
modes (k ~ 0.019 /i Mpc -1 ) from each band of multipoles for 
the highest waveband in Fig. 2. The scale of the histogram is 
marked on the right axis. For lower wavebands, the distribution 
moves to lower multipoles. The multipole number I is related to 
9 M by Z~3607<? M . 



any form of photometry errors), the rms residual of rela- 
tive flux error will be set by the intrinsic rms fluctuation 
of the projected galaxy number density due to cluster- 
ing, which is well below 1% on scales above several de- 
grees, where the contribution to Nk (k ~ 0.019 ft. Mpc -1 ) 
peaks. For smaller fcs (larger scales), the contribution 
comes from lower multipoles, and the rms fluctuation of 
the number counts decreases further. This means that 
the relative flux error over several square degrees can be 
controlled to better than 1% by two-dimensional galaxy 
counts from the LSST, which is sufficiently accurate for 
measuring the three-dimensional power spectrum on the 
largest scales. 

Further improvement on the relative flux error is possi- 
ble by combining galaxy counts with multi-band photom- 
etry (e.g. Babbedgc, Whitaker, & Morris 2005) and H I 
and CO surveys. Since the Poisson noise in the galaxy 
counts is an order of magnitude lower than that caused 
by galaxy clustering, one can also divide the galaxies 
into groups of similar properties and compare them in 
one field with those in another to better determine the 
differential extinction. 

5. APPLICATIONS 

Determination of the mass power spectrum has many 
scientific applications. They include searching for dark 
energy fluctuations, determining the dark energy equa- 
tion of state, determining the sum of neutrino masses, 
and probing the primordial power spectrum of fluctua- 
tions produced during inflation. Here we consider two ap- 
plications: determining the primordial power spectrum 
to probe inflation and determining the distance-redshift 
relation in order to probe dark energy. 

5.1. The Primordial Power Spectrum of Density 
Perturbations 



The mass power spectrum on large scales is a direct 
measure of the primordial fluctuations, and it provides 
a means to probe the generator of these fluctuations. 
Recent CMB observations have revealed some puzzling 
properties of the largest scales (Peiris et al. 2003; Efs- 
tathiou 2004; Schwarz et al. 2004; de Oliveira-Costa et al. 
2004; Eriksen et al. 2004b; Hansen et al. 2004; Land & 
Magueijo 2005; Jaffe et al. 2005). These peculiar features 
may very well have their origin in systematic error. How- 
ever, they have at the very least served to draw attention 
to the possibility that the very large-scale structure of 
the Universe may be more complex than in our simplest 
models. 

For specificity, let us focus on one of these puzzles and 
see how a large- volume photometric redshift survey could 
shed further light on the solution. In Peiris et al. (2003) 
it is noted that there are some unusual sharp features, 
or 'glitches', in the rise to the first peak of the WMAP 
temperature power spectrum. Peiris et al. (2003) also 
suggest that the glitches may be evidence of features in 
the primordial power spectrum. Adjusting two parame- 
ters of an inflation potential model with a sharp break, 
or 'step' (e.g. Gottlober, Miiller, Starobinsky 1991), they 
find a significantly better fit to the power spectrum than 
if they assume a featureless power law. The step inflation 
model leads to a new physical length scale (the size of the 
Horizon as the inflation scalar field crosses the feature in 
its effective potential) that imprints itself on the power 
spectrum as a series of peaks and troughs. 

Peiris et al. (2003) point out that the feature in the 
primordial power spectrum is detectable in future large- 
volume redshift surveys. One can indeed see in Fig. 2 
that the statistical power is there to map out the wiggles 
with a high signal-to-noise ratio. 

These features are much more pronounced in the mat- 
ter power spectrum than in the temperature power spec- 
trum. The CMB projection from three dimensions to two 
dimensions makes them less prominent. In general it is 
interesting to compare the relative abilities of the CMB 
and galaxy surveys to constrain the primordial power 
spectrum. 

Before turning to forecasts from CMB data though, 
let us first refine our forecasts from galaxy survey data 
at k > 0.01 ft, Mpc -1 . As seen in Fig. 4, one has to 
apply a large correction for photometric redshift errors 
and, to a much less extent, redshift distortions to recover 
the matter power spectrum at k > 0.01 ft Mpc -1 . An 
incorrect probability distribution for the redshift errors 
can introduce significant errors in the recovered matter 
power spectrum. Even if the correction is accurately de- 
termined, one will still be limited by the shot noise at 
k > 0.1 ft Mpc -1 , which is boosted by the correction 
along with the observed galaxy power spectrum. 

To avoid sensitivity to uncertain and large corrections 
for suppression of power due to photometric redshift er- 
rors, we can simply discard modes with a lot of suppres- 
sion in the radial direction. That is, as discussed in Sco 
& Eisenstein (2003) and Glazebrook k Blake (2005), we 
can discard modes with fc|| above some critical value. The 
approximate mode counting equation (4) becomes 
J k 2 AkV/2ir 2 fey < k* { 
Nk = { kk*AkV/2n 2 fc|| > k*, ( U > 

where a cubic geometry and the plane-parallel approxi- 
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Fig. 7. — Forecasts of sample variance errors on the primordial 
matter power spectrum. For the LSST, we set z max = 1-0 (dotted 
line) and 2.5 (solid line), and the mode counting follows equation 
(14) with k* = 0.01 hMpc" 1 . The forecast for the CMB (open 

squares) includes both temperature and polarization information, 
and it is taken from Hu & Okamoto (2004). Both forecasts assume 
a binning of Ak = 0.05A:. 



mation are assumed. One would choose fc^ to be small 
enough so that residual errors in the power spectrum can 
be tolerated after correcting for photometric rcdshift er- 
rors and redshift distortion. We set fcj| = 0.01 /iMpc -1 

and assume that the residual error is negligible for fc|| < 
fc* 

Note that at low fc the number of modes scales with k 
as a three-dimensional survey, increasing as fc 3 if Ak oc k, 
and that above fcjj the scaling is that of a two-dimensional 
survey. Thus, although the CMB will always be the best 
on length scales larger than those that fit in the galaxy 
survey, a galaxy survey can compete, at least in terms of 
the statistical weight, with the CMB on sufficiently small 
scales. If the cross over scale is large enough, galaxy bias 
modeling will be simple enough to exploit this statistical 
improvement. 

We present in Fig. 7 the sample variance errors of the 
primordial matter power spectrum for the LSST with 
•Zmax = 1-0 (dotted line) and 2.5 (solid line). Also in- 
cluded is the error forecast for the CMB (open squares) 
using both temperature and polarization information 
(Hu & Okamoto 2004). It is very encouraging that the 
LSST can almost match the best-case-scenario CMB re- 
sults with z max = 1.0 and even provide a better measure- 
ment of the primordial power spectrum with z max = 2.5. 
The cross over scale referred to in the previous paragraph 
is indeed at very large scales. For the z max = 2.5 case, it 
is at fc < 5 x 10" 3 /iMpc" 1 . 

The binning in Fig. 7 is for Afc = 0.05fc for ease of com- 
parison with Hu & Okamoto's results. However, the spec- 
tral resolution of the galaxy survey is 7.5 x 10" 4 h Mpc -1 
(•Zmax = 2.5) with exceptions for the first few modes of 
each multipole, which have slightly coarser resolutions. 
Thus for k < 0.015 /iMpc -1 the binning assumed in 
Fig. 7 is too fine. To avoid correlations, that would oth- 
erwise increase the errors above the values plotted in the 
figure, one could bin more coarsely. This coarser binning 
would reduce the level of error below what is plotted in 



the figure by a factor of -\/0.05fc/Afc, although the num- 
ber of points would, of course, be reduced also. 

Hu & Okamoto (2004) point out that the CMB power 
spectrum errors in Fig. 7 arc highly correlated and that 
there are certain linear combinations of the power spec- 
trum measurements with much smaller errors. Approxi- 
mately 50 principal-component modes of the covariancc 
matrix can be measured to percent level precision. The 
best-determined principal-component modes are combi- 
nations of Fourier powers at fc > 0.08 h Mpc" 1 . 

Our forecast for how well the primordial power spec- 
trum can be determined from galaxy survey data is done 
in the limit of no uncertainty in r(z), g(z), b(z), u) m 
and cjb- We expect that CMB data from Planck and 
the galaxy survey data itself can be used to reconstruct 
these well enough that the uncertainties do not quali- 
tatively change our error forecasts. The Hu & Okamoto 
(2004) forecast does include the effect of these uncertain- 
ties which are more important for the CMB because of 
the projection from three dimensions to two dimensions 
and the greater prominence of the acoustic oscillation 
features. 

We must mention two further caveats. First, for the 
lowest fc values the spherical geometry slightly increases 
the variance beyond what one gets for a cubic geome- 
try as shown in Fig. 2. Second, whether one can neglect 
errors in the correction for photometric rcdshift power 
suppression at fc < 0.01 ft. Mpc -1 remains to be demon- 
strated. 

5.2. Standard Rulers For Determining r(z) 

In the previous subsection we considered how well the 
primordial power spectrum could be determined with 
r(z), g(z), b(z), LUb and uj m known perfectly. Now we 
consider how well r(z) can be determined if we assume 
the primordial power spectrum to be featureless, as ex- 
pected in the simplest models of inflation. If the power 
spectrum remained featureless, it would be impossible 
to reconstruct r{z) from galaxy clustering data, but fea- 
tures imprinted by the evolution of this spectrum can be 
used as standard rulers for determination of the angular- 
diameter distance. 

In linear perturbation theory there are two length 
scales that become imprinted on the matter power spec- 
trum. The larger scale feature is the peak near fc = 
0.02/i Mpc" 1 that can be seen in Fig. 3. The location 
of this feature depends on the size of the horizon at the 
epoch of matter-radiation equality. Assuming a standard 
radiation content, this is set by the matter density, w m , 
which can be determined to sub-percent accuracy with 
Planck (Eisenstein, Hu, & Tegmark 1999). The smaller 
scale feature is a series of peaks and troughs in the matter 
power spectrum due to acoustic oscillations the baryons 
undergo in the pre-recombination plasma (Peebles & Yu 
1970; Bond & Efstathiou 1984; Holtzman 1989; Hu & 
Sugiyama 1996). The length scale here is the comoving 
size of the sound horizon at last-scattering. This de- 
pends on cu m and Ub which can both be determined well 
from CMB observations. Thus, both of these features 
can serve as CMB-calibratcd standard rulers, which one 
can use to determine the comoving angular-diameter dis- 
tance, r(z) (Eisenstein et al. 1998; Cooray et al. 2001). 

Since r(z) = c J dz' / 'H(z') these measurements can 
be used to map out the history of the expansion rate and 
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thereby provide constraints on the dark energy. 

First we will consider the baryonic oscillations. The 
lowest k peak in the series of baryonic oscillations has 
already been observed in spectroscopic redshift surveys 
(Eisenstein et al. 2005; Cole et al. 2005). Several papers 
have studied how well r(z) and cosmological parameters 
can be determined from future surveys, both spectro- 
scopic and photometric (Blake & Glazebrook 2003; Hu 
& Haiman 2003; Linder 2003; Seo & Eisenstein 2003; 
Blake & Bridle 2004; Dolney, Jain, & Takada 2004; An- 
gulo et al. 2005; Glazebrook & Blake 2005; Linder 2005; 
Zhan & Knox 2005). Prospects for controlling non- linear 
evolution and galaxy biasing look very good (White 2005; 
Seo & Eisenstein 2005). Photometric surveys over larger 
volumes can provide constraints that are competitive 
with spectroscopic surveys over smaller volumes. 

If one counts the modes using equation (14), the shot 
noise will continue to be negligible at k > 0.1/iMpc -1 
and the dominant errors there will come from the un- 
certainties of the correction for photometric redshift er- 
rors, nonlinear redshift distortion, nonlinear evolution, 
and galaxy bias. Despite the complexities, these errors 
do not produce oscillating features in the power spec- 
trum. Thus it is possible to use a photometric redshift 
survey to measure the angular-diameter distance accu- 
rately. 

We now turn to a discussion of the large-scale feature. 
Current galaxy surveys cannot be used to measure the 
power spectrum on scales larger than k = 0.02 /iMpc -1 
due to their limited survey volume. The LSST, on the 
other hand, will be able to probe scales with wavenum- 
bers as small as several thousandths /iMpc -1 . At the 
expense of increased sample variance error, one may 
measure the matter power spectrum in several redshift 
bins with the LSST. For example, Fig. 7 suggests that 
4 equal-volume bins from z = to 2.5 will enable us 
to measure the matter power spectrum to roughly 10% 
around k = 0.01 Mpc' 1 . The angular scale to which 
this feature projects depends on the angular-diameter 
distance to each redshift shell, enabling one to determine 
the angular-diameter distance. 

Generally, constraints from the large-scale feature will 
be weaker than from the baryonic oscillations for two rea- 
sons: first, the latter are much sharper and therefore less 
tolerant of the horizontal shifts induced by a re-scaling 
of distances and second, the sample variance errors are 
larger on larger scales. 

However, even this broad feature can lead to powerful 
distance determinations. Seo & Eisenstein (2003) find 
that the angular-diameter distance determinations using 
the large-scale feature have 2 to 3 times larger errors. 
Thus, though subdominant in their constraining power, 
the large scale feature offers a useful independent check 
on distance determination from the baryonic oscillations. 

With the LSST and CMB priors from the Planck mis- 
sion, one can achieve errors of ~ 1% on r(z), 0.10 
on wo, and 0.25 on iy a using both the baryon oscilla- 
tions and broadband feature (Zhan & Knox 2005), where 
the dark energy equation of state is parametrized as 
w(z) = w + Wg.[l - (1 + z)- 1 ]. 

6. DISCUSSION AND CONCLUSIONS 

For the reconstruction of the matter power spectrum 
P(k) to work, we must know the comoving distance r(z), 



the linear growth function g{z), the galaxy bias b(z), and 
the rms photometric redshift error a z adequately well. 

One may improve the knowledge of the galaxy bias 
on large scales by comparing galaxy statistics with mass 
statistics inferred from weak lensing (Hoekstra et al. 
2002; Seljak et al. 2005; for a joint analysis, see Hu & 
Jain 2004), by constructing realistic theoretical models 
(e.g. Peacock & Smith 2000), and by studying galaxies 
in large volume hydrodynamical simulations that include 
star formation and feedback (e.g. Weinberg et al. 2004). 
Available observational constraints on bias (e.g. Tegmark 
et al. 2004) constrain it on scales smaller than those of 
interest here, but since the LSST is designed for weak 
lensing, it can help to measure the galaxy bias on large 
scales for galaxies associated with the lensing mass field. 

The uncertainties in r(z) and g(z) can be mitigated 
by assuming a prior for the cosmological model, but at 
the same time we lose the ability to measure the matter 
power spectrum model-independently. However, since 
cosmological information is encoded in the matter power 
spectrum through multiple channels in addition to r{z) 
and g(z), cosmological parameters can be estimated from 
the reconstructed matter power spectrum as well. This 
will provide a self-consistency check for the choice of the 
parameters. 

Photometric redshift errors strongly suppress the raw 
galaxy power spectrum (\Di n \ 2 ). Thus, it is crucial to 
determine the correction for photometric redshift errors. 
We have adopted a form of this correction from Heavens 
& Taylor (1995), but other forms also exist (e.g. Ballingcr 
et al. 1996). They should be tested against A-body sim- 
ulations to improve the accuracy. Meanwhile, it appears 
possible from Fig. 4 that one can infer the correction by 
comparing the galaxy power spectrum of different multi- 
poles. One may also discard high fc|| modes to subdue the 
errors due to uncertainties in the correction (see Fig. 7) . 

The extinction is yet another challenge to measuring 
the matter power spectrum on very large scales. We have 
demonstrated that with 250,000 projected galaxies per 
square degree the LSST can actually determine the dif- 
ferential extinction to better than 1% from galaxy counts 
alone. This means that not only can we be sure that the 
signal is truly coming from clustering of galaxies but we 
can also combine the galaxy counts with multi-band pho- 
tometry (e.g. Babbedge et al. 2005) and H I and CO maps 
to produce a more accurate reddening map, which will 
be useful for other observations such as CMB surveys. 

Despite all the difficulties, we find that a deep and 
wide photometric redshift survey, such as can be done 
with the LSST, can measure the power spectrum of pri- 
mordial fluctuations on very large scales. The shot noise 
on these scales is so low (see Fig. 4), that it is possible 
and beneficial to divide the galaxies into different types 
to check for consistency and to indirectly test the scale- 
independence of the galaxy bias on large scales. 

The large-scale primordial power spectrum measured 
from galaxies will complement the CMB experiments and 
provide valuable insights on inflation. For example, it 
can be used to examine the possibility of a step infla- 
tion potential (Peiris et al. 2003). The large-scale peak 
and intermediate-scale baryon oscillations in the matter 
power spectrum will provide useful constraints on cosmo- 
logical parameters including the dark energy equation of 
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state parameters. They can also reduce the uncertain- 
ties of the photometric redshift bias (Zhan & Knox 2005), 
and thus prevent the redshift bias from severely degrad- 
ing the constraining power of the weak lensing tomog- 
raphy on the dark energy equation of state parameters 
(Huterer et al. 2005; Ma, Hu, & Huterer 2005). 
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APPENDIX 
SPHERICAL HARMONIC ANALYSIS 

For a survey that covers a large fraction of the sky, the spherical coordinate system is a more suitable basis to work 
with. The spherical harmonic analysis has been applied to the IRAS galaxy catalog (Scharf et al. 1992; Fisher, Scharf, 
& Lahav 1994), the PSCz galaxy catalog (Tadros et al. 1999), and, recently, the 2dFGRS catalog (Percival et al. 2004). 
We adopt the formulation in Heavens & Taylor (1995) and describe it here briefly. 

The spherical harmonic decomposition of a quantity ^4(r) is 

Ai mn =ci n [ w(r,fc hl )A(r)j i (fc i „r)y i ^(f)d 3 r, (Al) 
Jo 

where an extra weight function u>(r, fc;„) is inserted to minimize the variance in the power spectrum measurement. The 
discrete wavenumber ki n is the nth root of the boundary condition dji(fer)/dr| rmax = 0. The normalization q„ satisfies 
c hi = / 3? (ki n r)r 2 dr . One can relate the observed three-dimensional galaxy overdensity Di, nn to the underlying mass 
overdensity Si mn through 

Dlmn ^ ^ Sinn' Qln'n" &lmn" 3 0^-^) 
n'n" 

Sinn' —cinCin' / 7= ji{ki n r)ji{ki n ' y)rdr ydy , (A3) 



2na z 

Qinn' =cinCin< J n s (r)g(z)^b(z)w(r, k)ji{kin'r)ji{ki n r) + 

O - 6 M d d -1 

-J*^JL[ w (r, k)Mki nr )]^Mk M r)}^dr, (A4) 

where a z is the photometric redshift error in comoving distance units, Q m (z) is the ratio between the cosmic matter 
density to the critical density at redshift z, and we have applied the approximation dlng/dlna ~ fi"; 6 (Lahav et al. 
1991). For simplicity, we assume a full-sky survey with no angular variation in the selection function, so that the 
galaxy distribution n g (r) = n g (r)[b(z)g(z)S(r) + 1], in which S(r) is scaled to z = 0. The term Si nn ' would account for 
the nonlinear redshift distortion, i.e. the pair-wise velocity dispersion, if there were no error in redshift measurements. 
Since the photometric redshift error is far greater than the pair-wise velocity dispersion, the latter is neglected. The 
term Qi nn ' is due to the linear redshift distortion (see also Scoccimarro 2004). With (8imnSi' m >n>) ~ ^(^nj^'t'C'i 
one can show that 

(\D ln \ 2 ) = (\D lmn \ 2 ) = ]T U? nn ,P(ktn>) + N ln , (A5) 

n' 

where the window function U? , and shot noise Ni n , respectively, are given by 

Uinn' = ^2 S lrm"Qln"n' and N ln = cf n n g (r)w 2 (r, k)jf {k ln r)r 2 dr. (A6) 

n" 

The weight function w(r, k) is introduced in equation (Al) to reduce the variance in the power spectrum estimator 
(Fcldman et al. 1994; Heavens & Taylor 1995). Designing an optimal weight function is beyond the scope of this 
work. Here, we only demonstrate the performance of the estimator with a weight function w(r,k) — P(k)/[1 + 
n g (r)g(z)b(z)P(k)] (see also Feldman et al. 1994; Heavens & Taylor 1995). Fig. 5 shows several examples of the 
window function U 2 nn , . The fiducial cosmological parameters are given by the WMAP results (Spergel et al. 2003), 
and the matter power spectrum is calculated using the fitting formula from Eisenstein & Hu (1999). We model 
the evolution of an overall bias with b(z) = 0.84z + 1 (e.g. Weinberg et al. 2004). The symbols mark the discrete 
wavenumbers of the radial modes. The window function U 2 nn , is fairly narrow for k < 0.02 ftMpc -1 , except for Uf ln , 
(not shown), which are discarded in the mode counting. The results of other multipolcs share the same characteristics, 
except that for higher multipoles the modes move toward large wavenumbers and vice versa. 
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